Dragging two-dimensional discrete solitons by moving linear defects 
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We study the mobility of small-amplitude solitons attached to moving defects which drag the 
solitons across a two-dimensional (2D) discrete nonlinear-Schrodinger (DNLS) lattice. Findings 
are compared to the situation when a free small-amplitude 2D discrete soliton is kicked in the 
uniform lattice. In agreement with previously known results, after a period of transient motion the 
free soliton transforms into a localized mode pinned by the Peierls-Nabarro potential, irrespective 
of the initial velocity. However, the soliton attached to the moving defect can be dragged over 
an indefinitely long distance (including routes with abrupt turns and circular trajectories) virtually 
without losses, provided that the dragging velocity is smaller than a certain critical value. Collisions 
between solitons dragged by two defects in opposite directions are studied too. If the velocity is small 
enough, the collision leads to a spontaneous symmetry breaking, featuring fusion of two solitons into 
a single one, which remains attached to either of the two defects. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.75.-b 



I. INTRODUCTION 

The discrete nonlinear Schrodinger (DNLS) equations 
constitute a vast class of systems which are profoundly 
interesting in their own right, and serve as important 
physical models for nonlinear optics lj, matter waves in 
Bose-Einstein condensates (BECs) J2J], and in other con- 
texts [3[. In particular, soliton solutions to the DNLS 
equation in one, two, and three dimensions (ID, 2D, and 
3D) represent fundamental localized modes in discrete 
media. Experimentally, ID and 2D quasi-discrete soli- 
tons have been created in nonlinear optical systems of 
several types, see a comprehensive review in Ref. [1|. 

Local defects are important ingredients of DNLS mod- 
els. They are interesting as additional dynamical ele- 
ments of the lattices [3], and have direct physical re- 
alizations. In particular, they may describe various 
strongly localized structures in photonic crystals Q , such 
as nanocavities [5[, micro-resonators [6], and quantum 
dots [7]. In the context of BEC, local defects can be eas- 
ily created and used for manipulating the condensate by 
means of focused laser beams. This technique has made 
it possible to create optical tweezers for trapping and con- 
trollable transfer of condensates by moving optical traps 
j8| . A scheme for nonlinear optical tweezers that can ex- 
tract solitons from a linear reservoir was proposed too 
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In a more general context, controllable transport of 



solitons in periodic and quasiperiodic media, including 
the ultimate case of discrete lattices, play s an important 
role in various applications, see Refs. [lOj and references 
therein. Many attempts have been made to devise a sim- 
ple way of transferring solitons from one position to an- 
other with minimal losses. In particular, dragging gap 
solitons in BEC, embedded into an optical lattice, by a 
moving defect to which the soliton is attached, was pro- 
posed as a means of the transport in Ref. [ll[ . Recently, 
manipulations of a BEC vortex by a localized impurity 
representing a focused laser beam was considered in Ref. 

m. 
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The subject of the present work is the dragging of 
discrete solitons in 2D lattices by an attractive defect 
which plays the role of the mover. The defect is shaped 
as a Gaussian of a finite width. After introducing the 
model and recapitulating some relevant results for static 
trapped modes in Sect. 2, we consider moving solitons in 
Sect. 3. First, the problem of the immobility of free 2D 
DNLS solitons is briefly revisited, and then new results 
are reported for the transfer of solitons by the moving at- 
tractive defect. Both simple rectilinear dragging routes, 
and more sophisticated ones, in the form of square, rhom- 
bic, and circular closed trajectories, are considered. It 
is concluded that the trapped solitons may survive the 
dragging over an indefinitely long distance, if the drag- 
ging velocity does not exceed a certain critical value. The 
paper is concluded by Sect. 4. 



II. THE MODEL AND STATIONARY 
SOLUTIONS 



B. Existence curves for stationary solitons and 
stability 



A. The formulation 

We consider the following model based on the 2D 
DNLS equation with a local defect: 



-JA2Un,m-\-V ni m(t)Un,m-\-^\Un,r. 



o, (1) 



where the overdot stands for the time derivative, 

^2^n, m — ^n,m+l ~r ^n,m— 1 ~r ^n+l,m ~r ^n— l,m ^^n^m IS 

the 2D discrete Laplacian, the coupling constant of the 
lattice will be fixed by scaling, J = 1, and a > is the 
coefficient of the on-site self-attraction. Discrete function 
V^ m account for the linear defect, taken in the form as 
a Gaussians profile, 

V ni m(t) = vexp {- [(n - n(t)f + (m - m(t)) 2 ] /A 2 } , 

(2) 
with respective strength v, width A, and coordinates of 
the (moving) center, n(t) and m(t). In this notation, pos- 
itive and negative strengths correspond to the attractive 
and repulsive defect, respectively. 

Looking for stationary solutions in the case of the qui- 
escent defect, with m(t) = n(t) = in Eq. (J2j), 

U n ,m = U n , m exp(-iuot), (3) 

we arrive at the nonlinear eigenvalue problem, 

uU n ,m + A 2 Z7 n , m + V nirn U n ,m + 0"|^n, m\ U n , m = 0, (4) 

for real frequency uo and the profile of the stationary dis- 
crete mode, Um,n^ which may be complex, in the general 
case. In the absence of the defect, the linearized version 
of Eq. (J2J) gives rise to the dispersion relation for linear 
modes Um,n ~ ex P (ikm + iqn), 



4 [sin 2 (fc/2) + sin 2 (q/2)] , 



(5) 



which features the phonon band, < uo < 8. Above 
and below the band, full nonlinear equation ^ gives 
rise to nonlinear modes, described by respective curves 



N(uo), with N = J2 Ti 



being the norm (power) 



of the mode. The important difference between 2D and 
ID settings is that, in the latter case, the fundamental 
single-peak mode (alias the Sievers-Takeno discrete soli- 
ton) admits the limit of TV — )> 0, while all the 2D solitons 
are bounded by a critical norm, N cr , below which they 
do not exist [3]. Accordingly, an abrupt derealization 
(decay) of discrete 2D solitons was predicted in the case 
when the inter-site coupling constant exceeds a certain 
critical value [13J |. It is worthy to mention that, intro- 
ducing a defect with the linear and nonlinear components 
and varying their strengths, or the lattice coupling con- 
stant, in time (in the spirit of the "soliton management" 
techniques |l3|). one can change the critical norm, 7V C r , 
and thus control the transition to the derealization [15j . 



Curves N(uo) for the discrete solitons obtained from nu- 
merical solutions of Eq. (J2J), based on the continuation 
from the anti-continuum limit [the one corresponding to 
J = in Eq. (pQ)], using the standard Newton's itera- 
tion routine, are displayed in Fig. [TJ in the presence and 
in the absence of the defect. In panel QJa) we compare 
three curves: the black one for the case of a = 1 without 
the defect, and the blue and red curves for a = 1 and 
a = 0.25 with the defect of amplitude v = 1 and width 
A = 1 [the red curve is obtained from the blue one by 
rescaling of the norm, in order to get a curve crossing 
the black one in a vicinity of uo — >• 0, see Fig. DJb)]. 
In the defect-free lattice, a family of stable solitons (in 
agreement with the Vakhitov-Kolokolov (VK) criterion, 
it obeys condition dN/duo < [3]) can be found in the re- 
gion of — oo < uo < uo cr w — 1. In this region, the discrete 
soliton is tightly localized, being immobile, due to the 
strong pinning to the underlying Peierls-Nabarro poten- 
tial |3| . The Peierls-Nabarro potential becomes weak for 
broad small-amplitude discrete solitons corresponding to 
uo — » 0. Note that the introduction of the positive defect 
can significantly reduce or completely suppress the insta- 
bility region, as shown by the blue and red curves in Fig. 
QJa) where only a small segment of the existence curve 
in the interval of —0.58 < uo < —0.52 remains unstable. 
Note also that the attractive defect removes the lower 
existence bound for the 2D discrete solitons, 7V cr , as the 
soliton goes over into a linear defect mode in the limit of 
7V^0. 
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FIG. 1: (Color online) (a) The existence curves for on-site 
symmetric solitons in the defect-free lattice (v = 0, the black 
line) , and in the presence of the quiescent defect © with v = 1 
and A = 1, for different values of the nonlinearity coefficient: 
a — 1 and a = 0.25 (the blue and red lines, respectively). 
The dashed part of the black curve corresponds to unstable 
solitons in the uniform lattice. Panel (b) shows a zoom of 
the dashed box from (a). Points A and B correspond to the 
solitons with frequencies to a — —0.1 and ujb — —0.17, and 
equal norms, Na — Nb- 

According to the previous analysis, at small \uo\ initially 
unstable discrete solitons in the uniform lattice may be 
stabilized by switching the attractive defect on [3|, fl5j | . 
To keep the norm of the soliton constant in this case, 
simultaneously with introducing the defect we decrease 




FIG. 2: Profiles of unstable and stable discrete solitons, which 
correspond, respectively, to points A and B of the existence 
curves in Fig. [TJb). 



the strength of the nonlinearity. As a result, we obtain 
the existence curve shown in red in Fig. Ufa), which 
features intersection with the original black line in the 
region of \cj\ <C 1. The stabilization corresponds to the 
transition from point A to B in Fig. [ljb), which follows 
the increase of the strength of the defect from v = to 
v = 1 and reduction of the strength of the nonlinear coef- 
ficient from a = 1 to a = 0.25. We have checked numeri- 
cally that the transition between these two configurations 
could be performed even instantaneously in time, leading 
to a very small radiation loss. Examples of the profiles of 
the discrete solitons found without and with the defect 
are displayed in Fig. [2j Both solutions have the same 
norm, while solitons A and B belong, respectively, to the 
unstable and stable branches of the existence curves, cf. 

Fig.m 



III. MOVING SOLITONS 




FIG. 3: Snapshots of the discrete soliton initially kicked in the 
horizontal direction (c m = 0), are shown at different times. 
The initial profile is taken from Fig. Efa). Parameters of 
the kick and the coordinate of the final pinned site, rif, are 
c n = 0.1, n/ = 7 in (a); c n = 0.2, n/ = 9 in (b), and c n = 0.5, 
n/ = 12 in (c). 



A. Dynamics of discrete solitons in the uniform 
lattice 



Proceeding to the consideration of traveling solitons, 
it is first relevant to recapitulate basic results concern- 
ing the ballistic motion of kicked solitons in the uniform 
lattice. We kick the discrete soliton by taking the ini- 
tial condition as Z7 n?m = (7^4 exp [— i(c n n + c m ra)/2], 

where Un,L is stationary configuration ([3]), and vector 
(c n , Cm) determines the strength and direction of the ini- 
tial impulse. Figure [3] demonstrates that the broad small- 
amplitude soliton from Fig. [2fa), if kicked in the n di- 
rection by the impulses of different strengths, c n = 0.1, 
0.2, and 0.5, is transformed into a tightly localized peak, 
which keeps nearly the entire initial norm and comes to 
the halt at a site with coordinate n/. Note also that 
the initial broad soliton was unstable, while the final tall 
peak (with almost the same norm) represents a stable 
soliton, as per the black curve in Fig. [TJ 

This and other simulations corroborate the known fact 
that, under the action of the kick applied in any direc- 
tion relative to the underlying 2D lattice, the broad 2D 
DNLS soliton starts to move, but eventually stops after 
having traveled a finite distance (on the contrary to the 
ID DNLS model, where the kicked soliton may travel in- 
definitely long [3, [l6| ) . The harder the initial shove c n , 
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FIG. 4: (Color online) The evolution of the peak density for 
three different values of the kick: c n = 0.1 (solid black), 
c n — 0.2 (dashed red), and c n = 0.5 (dotted blue). Other 
parameters are as in Fig. [3] 



the longer distance n/ is passed by the discrete soliton 
before the stoppage (cf. Figj3]). However, an excessively 
strong kick induces strong perturbations in the shape of 
the soliton, under the action of which it starts to radi- 
ate, leading to a significant loss of the norm. In Fig. 2J 
the peak density of the moving soliton, determined as 
Umax = max(|7i n>m | 2 ), is displayed for different strengths 
of the initial kick. It is seen that, after coming to the 
halt, the initially broad soliton transforms into a highly 
localized mode whose amplitude features gradually fad- 
ing oscillations. Similar results were obtained for the 



propagation of the discrete soliton initially kicked in the 
diagonal direction, with c n = c m . 

Thus, persistent motion of solitons in the uniform 2D 
DNLS lattice is impossible. This may be explained by the 
fact that, in the continual limit, the cubic self- focusing 
leads to the collapse of solitons; in the discrete setting, 
this implies the formation of tightly localized tall peaks, 
which are strongly pinned to the lattice, being therefore 
immobile. On the other hand, it is known that motile 2D 
solitons are possible in discrete media with weaker non- 
linearities, which do not lead to the collapse in the con- 
tinual limit, viz., saturable [17] and quadratic [18| on-site 
nonlinear terms. In the 2D lattice with the combination 
of self-focusing cubic and defocusing quintic terms (this 
combination of the nonlinearities does not give rise to 
the collapse in the continuum medium either), the kicked 
soliton may perform a long run, but eventually it comes 
to a halt, because of radiation losses (l9l. 



B. Dragging the discrete solitons by the moving 
defect along the linear trajectory 

Forced motion of solitons attached to moving defects 
is the central topic of this work. To demonstrate a typi- 
cal realization of this scenario, we take the initial soliton 
from Fig. Efb), and drive the defect in the lattice plane 
along a rectilinear rout: n(t) = n s +c n t, m(t) = m s +c m t, 
where c n and c m are the velocities in the n and m di- 
rections. A set of typical examples is displayed in Fig. 
[5j for dragging the discrete soliton in the n direction 
with different velocities, c n = 0.1, and 0.5. This fig- 
ure demonstrates that the soliton can be transferred, 
virtually without any loss of the norm, over indefinitely 
long distances, if the velocity, which is applied instan- 
taneously, is small enough, allowing the soliton to per- 
manently adjust itself, in an adiabatic manner, to the 
positions passed in the course of the motion [Fig. [5ja)]. 
In more quantitative terms, this may be explained by 
the comparison of the dragging velocity with the vecto- 
rial group velocity produced by dispersion relation (J5j): 
v = duj/d{k,q} = 2 {sin fc, sin q}. If the dragging ve- 
locity is not small enough in comparison with the group 
velocity, the adiabatic adjustment of the moving soliton 
is not possible, and the soliton will be destroyed by the 
emission of radiation waves escaping at the group veloc- 
ity. 

To further illustrate peculiarities of the dynamics of 
the driven soliton, its peak density is displayed in Fig. 
[6j as a function of the rescaled time, for different drag- 
ging velocities. The oscillations correspond to transfor- 
mations between on-site (maximum peak density) and 
off-site (minimum peak density) configurations of the dis- 
crete soliton. At small dragging velocities, the total norm 
remains practically constant after a small initial loss of 
the norm due to the radiation, while the large velocity 
(c n = 0.5 in Fig. [6]) causes a significant initial loss. In the 
latter case, the subsequent dynamics shows conservation 



of the total norm. Similar results (not displayed here) 
have been obtained for the driven motion of the soliton 
along the lattice diagonal. We do not aim here to exactly 
identify the critical velocity for the transition from the 
movable solitons to immovable ones, as the transition is 
somewhat fuzzy, going through an intermediate region 
where the radiation loss suddenly starts to increase with 
the growth of the driving velocity. 




FIG. 5: Snapshots of the discrete soliton, dragged by the 
moving defect, at different times. The initial profile, as well 
as parameters of the defect, are the same as in Fig. [2jb) 
(v = A = 1). The velocity of the defect is c n — 0.1 in (a), 
and c n = 0.5 in (b). The other component of the velocity is 
absent, c m = 0. 




FIG. 6: (Color online) The evolution of the peak density for 
two different values of the dragging velocity: c n = 0.1 (solid 
black), c n = 0.5 (dashed red). Other parameters are the same 
as in Fig. \5\ 



C. Dragging the soliton along closed trajectories 

To demonstrate more sophisticated kinetic effects, in 
Fig. [71 we display dragging the soliton along a closed 
trajectory formed by four segments aligned with bonds 
of the lattice. Typical results are also displayed in Fig. 
[8] for the driven motion along a closed route formed by 
four diagonal segments. Figures [3 and [8] suggest that, 



in the case of the abrupt changes in the direction of the 
driven motion, the velocity of the dragging defect is a 
crucial parameter, leading to strong losses at large veloc- 
ities. Indeed, the comparison of Figs. EJb) Efc) to Fig. 
0D(b) demonstrates that relatively high velocities produce 
a much more destructive effect on the solitons dragged 
along the closed trajectories than on their linearly driven 
counterparts. On the other hand, for sufficiently small 
velocities [c n = 0.1 in Figs. [TJ^a) and |5fa)], the soli- 
ton can be dragged along complex routes in the virtually 
intact form. Continuing the study in this direction, we 
have also considered driven circular motion of the soliton, 
as shown in Fig. [9j where the discrete soliton is safely 
dragged along the ring trajectory with radius R = 15 and 
angular velocity 0.01/7T. 
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FIG. 7: The same as in Fig. \5\ but for dragging the soli- 
ton along the square formed by the lattice axes. Here the 
corresponding velocities along the n and m directions are 
\c n \ = \cm\ = 0.1 in (a) and \c n \ = \c m \ = 0.5 in (b). In 
panel (c), the time evolution of the peak density Umax for 
both cases is shown. 




FIG. 8: The same as in Figs. [TJa) and (b), for dragging the 
soliton along the rhombus formed by diagonals of the lattice. 
Here \c n \ = \c m \ = 0.1, 0.3, and 0.5, in (a), (b) and (c), 
respectively. 
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D. Collisions between discrete solitons dragged by 
defects 



Another natural problem is collisions between defect- 
dragged solitons. First, we consider the head-on collision 
between identical solitons dragged toward each other. 
Colliding, they temporarily merge into a "lump" with the 
double norm, roughly adjusted to the defect with doubled 
strength. Continuing the simulations and monitoring at- 
tempts of the subsequent evolution of the merged lump, 
we observed two different scenarios, depending on the ve- 
locity of the dragging. If the velocity is higher then some 
critical value, after the separation of the two defects the 
lump splits back into two solitons, which are tantamount 
to those existing before the collision, see Fig. [TUJ 

The collision dynamics drastically changes if the veloc- 
ity is smaller than the critical value, see Fig. [TTJ Just 
before the collision, but when the two defects are still 



FIG. 9: Dynamics of the discrete soliton with initial profile 
taken from Fig. EJa) and initially placed at position (n, m) = 
(R, 0) with R = 15. The defect moves along the circular 
trajectory: rid = i?cos(0.01£/7r) and rrid = i?sin(0.01t/7r). 
The snapshots are shown with interval t — 250. In panel (b), 
the evolution of the peak density is shown. 



well separated in space (in the present case, at t « 750), 
one observes a symmetry breaking in the distribution of 
the local power (density) in the merging lump, which 
becomes unstable against oscillations of the total norm 
between the two defects. At the time of the collision, 
almost the same behavior occurs as in the previous ex- 
ample pertaining to the higher velocity. However, after 
the separation of the defects, the lump does not split, 
staying attached, as a single soliton, to one defect and 
moving with it, while the other defect remains "bare", as 
seen in Fig. [TTT c). 
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FIG. 10: The collision between two discrete solitons dragged 
by defects with velocities \c n \ = 0.05. Panels (a), (b) and 

(c) display profiles of the solitons in the initial state, at the 
collision moment, and in the final state, respectively. In panel 

(d) the respective evolution of the peak density is presented. 
The initial solitons are taken from Fig. [TJb). 
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FIG. 11: The same as in FigUHl but for velocities \c n \ = 0.01. 
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The symmetry breaking in the head-on collisions of 
the dragged solitons can be explained by the consid- 
eration of a linear-stability problem. To this end, we 
calculated the existence curves for static symmetric and 
asymmetric modes trapped by a pair of defects separated 
by different constant distances, and checked the linear 
stability of those modes, taking perturbed solutions as 
u n ,m = exp(-iut)[U n ,m + a n , m exp(At) + 6 n>m exp(A*£)], 



where A is an eigenvalue, U n , m is the unperturbed sta- 
tionary solution, and a n5m ,6 n?m are small perturbations. 
The so found unstable symmetric (dashed lines) and 
stable asymmetric (solid lines) families of the solutions 
are presented in Fig. [T2la), with the distance between 
the defects being S = 8 (black) and 10 (red) sites (the 
corresponding existence curves are practically overlap- 
ping). In Fig. fT2f b) dependences of the most unsta- 



ble eigenvalues (the real part of A) on frequency uo for 
the symmetric configurations with different distances be- 
tween the defects from Fig. fT2Ta) are presented. One 
can conclude that the symmetric configuration is unsta- 
ble against small perturbations (with the most unstable 
eigenvalue substantially diminishing with the increase of 
the separation between the defects). Although this in- 
stability pertains to the symmetric lumps supported by 
the pair of immobile defects, it provides an explanation 
to the instability of the merged symmetric lump also in 
the case of the slowly colliding defects. 
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FIG. 12: (Color online) (a) Existence curves for the families of 
symmetric (upper lines) and asymmetric (lower lines) discrete 
solitons supported by a symmetric pair of the defects. Solid 
and dashed portions of the curves correspond to stable and 
unstable solutions. Black and red lines (mostly overlapping) 
pertain to the configurations with the distance between the 
defects fixed to be 5 = 8 and 10 sites, respectively, (b) The 
frequency dependence of the real part of the two most un- 
stable eigenvalues for the symmetric solutions from (a). The 
black (+) and red (■) symbols correspond to the distances 
between the defects 5 = 8 and 10, respectively. 

As mentioned above, the symmetry breaking occurs 
between the colliding solitons driven by the slowly mov- 
ing defects. This fact can be explained by comparing the 
time necessary for the development of the instability of 
the lump produced by the collision and the collision time. 
The time of the mutual passage of the solitons attached to 
rapidly moving defects is smaller than the time required 
for the onset of the instability. In this case, the colli- 
sion is quasi-elastic, and it does not cause the breaking 
of the symmetry between the solitons, while for smaller 
velocities the collision time is sufficient for the develop- 
ment of the instability and ensuing symmetry breaking, 
which may eventually transform the two solitons into one, 



pinned to a single defect. Similar arguments were pre- 
viously used to explain the transition from quasi-elastic 
collisions to symmetry-breaking ones in several models 
supporting ID solitons in continual models [20j . 

We have also considered the collision of solitons with 
a mismatch in direction m perpendicular to the veloci- 
ties (the "aiming parameter"), £ = \m\ — rri2\ , where re- 
defines the initial coordinate of the j-th soliton in the di- 
rection of m. In particular, for the same small velocities 
as in Fig. [TTJ with two different values of the mismatch, 
£ = 8 and £ = 10, the results are shown in Fig. [13j 
It is seen that the character of the collision can be con- 
trolled by varying the mismatch, cf. Ref. [21(, where 
the collisions were simulated between 2D solitons in the 
framework of a continual dissipative equation. 



IV. CONCLUSION 

In this work, we have elaborated a scenario for the 
controlled transfer of 2D discrete solitons through the 
DNLS lattice. The main obstacle to the transport is the 
fact that 2D solitons in the uniform DNLS system cannot 
move persistently under the action of an initial kick, be- 
ing always braked by the underlying Peierls-Nabarro po- 
tential. Nevertheless, we have demonstrated that a rela- 
tively broad soliton can be dragged over indefinitely large 
distances, with virtually zero loss, by the moving attrac- 
tive defect, provided that the dragging velocity does not 
exceed a critical value. This is qualitatively explained by 
the comparison of the dragging velocity with the group 
velocity of the linear waves propagating in the uniform 
lattice. The stable driven motion is possible in any direc- 
tion, as well as along complex routes with corners, and 
along circular trajectories; however, the critical velocity 
is lower in the latter cases. Collisions between solitons 
dragged by two solitons in opposite directions were con- 
sidered too, with the conclusion that the two solitons 
spontaneously merge into a single one, which stays at- 
tached to either moving defect, if the collision velocity is 
small enough. 

These scenarios can be implemented (and used for var- 
ious applications) in the following form: a quiescent soli- 
ton may be prepared in the uniform 2D lattice, then the 
local attractive defect (s) may be induced [for instance, by 
laser beam(s) illuminating the corresponding BEC], effec- 
tively converting the free solitons into defect modes; next, 
the mode(s) may be transferred to a new position, as de- 
scribed above, and, eventually, the laser beam(s) may be 
switched off. Eventually, free solitons may be transferred 
according to this protocol, in the medium where these 
solitons are not motile by themselves. 

The analysis presented above may be naturally ex- 
tended in various directions, including the transfer of 
vortex solitons [22j (in the latter case, the trapping de- 
fect must be wide enough, for a sufficient overlap with the 
vortex; in fact, the defect itself may have a vortical struc- 
ture). A challenging generalization would be to develop a 



similar scenario for discrete solitons in 3D lattices, which may also be realized in terms of BEC. 
I 




FIG. 13: Color online) The collision with finite mismatch £. 
The snapshots in (a)-(c) and (d)-(f) correspond to £ = 8 and 
£ = 10, respectively. In panel (g), the respective evolution of 
the peak density is shown for £ = 8 (solid black) and £ = 10 
(dashed red). 
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